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ABSTRACT 

Aims. We seek to detect and quantify the effects of surface convection (granulation) on the line spectra of K-dwarfs as a first step 
towards a rigorous testing of hydrodynamic models for their atmospheres. 

Methods. Very high resolution (R 160,000 - 210,000), high signal-to-noise ratio (S/N > 300) spectra of nine bright K-dwarfs 
were obtained with the 2dcoude spectrograph on the 2.7 m Telescope at McDonald Observatory to determine wavelength shifts and 
asymmetries of Fe i lines. Spectra of the same stars acquired with the High Resolution Spectrograph (R ^ 120,000) on the 9.2 m 
Hobby Eberly Telescope were used as radial velocity templates to calibrate the wavelength scale of the 2dcoude spectra. 
Results. The observed shapes and positions of Fe i lines reveal asymmetries and wavelength shifts that indicate the presence of 
granulation. In particular, line bisectors show characteristic C-shapes while line core wavelengths are blueshifted by an amount that 
increases with decreasing equivalent width (EW). On average, Fe i line bisectors have a span that ranges from nearly for the weakest 
lines (residual core flux > 0.7) to about 75 m s for the strongest lines (residual core flux ^ 0.3) while wavelength shifts range from 
about -150 m s _1 in the weakest (EW a 10mA) lines to in the strongest (EW > 100mA) features. A more detailed inspection 
of the bisectors and wavelength shifts reveals star-to-star differences that are likely associated with differences in stellar parameters, 
projected rotational velocity, and stellar activity. While the former two are well understood and confirmed by our data, the relation 
to stellar activity, which is based on our finding that the largest departures from the expected behavior are seen in the most active 
stars, requires further investigation. For the inactive, slow projected rotational velocity stars, we detect, unequivocally, a plateau in the 
line-shifts at large EW values (EW > 100 mA), a behavior that had been identified before only in the solar spectrum. The detection of 
this plateau allows us to determine the zero point of the convective blueshifts, which is useful to determine absolute radial velocities. 
Thus, we are able to measure such velocities with a mean uncertainty of about 60 m s -1 for four of our sample stars. 

Key words, stars: atmospheres - stars: late-type - techniques: spectroscopic - sun: granulation 



1. Introduction 

The atmospheric structures of cool stars with convective en- 
velopes are inhomogeneous due to the interplay between the ra- 
diation field and the thermodynamic properties of the gas. Near 
the surface, radiation losses and the strong temperature sensitiv- 
ity of the continuum opacity induce a very steep temperature gra- 
dient, driving convection. The ascending cells of gas (granules) 
cool down quickly and the gas falls back into the stellar interior 
through filamentary structures (intergranular lanes). The down- 
flows are filamentary due to mass conservation and they pro- 
duce the buoyancy work necessary to support the upflows (e.g., 
IStein & Nordlundl [T998). This phenomenon, granulation, is di- 
rectly observed in spatially resolved images and spectrograms of 
the solar disk. For distant stars, and also for the Sun, the signa- 
tures of granulation can be detected in their sp ectra using precise 
measurements of absorption line profiles (e.g.. lGravll9 82. 2005; 
Gray & Tonet|l985| ; [Grav & Nage]|[l989[lDrayinsl ll987atT2008; 
Aspl und et all200ot (Allende Prieto et alj|2002lh 



The high degree of correlation between the intensity and 
velocity fields in inhomogeneous stellar atmospheres produces 
asymmetric lines (in particular when the stellar disk is unre- 
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solved), with characteristic shapes for the line bisectors (the 
midpoints of the horizontal segments across the wings of the 
lines). In addition, the wavelengths of spatially unresolved lines 
are blueshifted, owing to the larger contribution to the emergent 
flux from the brighter, hotter ascending gas. These Doppler shifts 
are often referred to as convective blueshifts. Note, however, that 
not all lines seen in a stellar spectrum suffer from these effects 
but only those that are formed in layers that are deep enough into 
the stellar atmosphere, where the correlation between tempera- 
ture and velocity fields due to granulation is strong. For example, 
molecular features, which form in cool high photospheric layers, 
are not expected to be blueshifted. 

Hydrodynamic simulations including the effect of the 
radiation field are necessary to better understand granula- 
tion. Two-dimensional, as well as three-dimensional radiative- 
hydrodynamic models have been compute d for the Sun (e.g., 
iNordlundl [19821; ISteffen & Frevtad fl 99 lh IStein & Nordlundl 
119981; [ Robinson et al.ll2003b IVogler||2004l) as well a s other stars 
(e.g. iNordlund & Dravinsi 11990b lAsplundl 120051; ICollet et al.l 
2006, 2007; TrampedacrT l2007l) . Despite the fact that, due to 
computational limitations, these models describe only a small, 
hopefully representative, piece of the stellar envelope and have a 
limited spatial resolution, they have been remarkably success- 
ful at reproducing t h e observed proper t ies of t he solar (e.g.. 



Drav ins et al.l I198U ISteffen & Frevtad 1 199 it lAsplund et al. 
2000) and stellar granulation (e.g.. iDravins & Nordlundl 1 1990; 
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Allende P rieto et alj|2002l) . These theory vs. observation com- 
parisons, however, have been performed for only a few stars and 
often using a limited set of spectral features. Thus, there is an ur- 
gent need for quantitative testing of the fundamental predictions 
of these hydrodynamic simulations in order to better understand 
stellar granulation and its impact on the observed spectra. 

Of particular interest is the impact of granulation on the de- 
termination of fundamental parameters and chemical composi- 
tions of stars. With few exceptions, most of these studies use 
"classical" one-dimensional model atmospheres, which are cal- 
culated adopting plane parallel geometry, hydrostatic equilib- 
rium, energy flux conservation, and local thermodynamic equi- 
librium (LTE), in addition to ignoring the effect s of magnetic 
fields , rotation, and chromospheric a ctivity (e.g.. lKuruczlll979l 
119931: lGustafssonetai]|1975i 120081) . Although state-of-the-art 
hydrodynamic simulations share with classical models many of 
these basic assumptions, they take into account the atmospheric 
inhomogeneities associated with the granulation phenomenon. 

Spectroscopic studies using classical model atmospheres 
have provided a wealth of useful information for general astro- 
physi cs and cosmology (se e, e.g., the reviews bv lWheeler et al.l 
1989 and lMc William! 1 997b and it is, therefore, very important to 
check whether the abundances determined using classical meth- 
ods are self-consistent. Exploring the impact of granulation (also 
referred to as "3D effects") on abundances and fundamental pa- 
rameters derived using classical methods, in particular for K- 
dwarfs (see below), is therefore crucial to improve our knowl- 
edge of cosmochemistry. 

Since the lifetimes of K-dwarfs (~ 20 Gyr) are larger than 
the age of the Galaxy (13.73±0.12Gvr, lHinshaw et alj2009b . no 
K-dwarf (and, evidently, no other cooler dwarf star) has evolved 
off the main sequence yet. Thus, as opposed to F and early G 
main-sequence stars, samples of K-dwarfs are free from biases 
in chemical evolution studies due to stellar death. Although M- 
dwarf and cooler star samples are also free from these biases, 
the analysis of their spectra is complicated by the abundance 
of molecules and condensates (dust, clouds), uncertain opac- 
ities and molecular equilibrium calculations, which make the 
model atmosphere comp utation difficult (e.g., [Hauschild t et al.l 
1 19991: iBurrows et al.ll2006h and the spectroscopic analysis ardu- 
ous (ill-defined continuum, atomic lines blended with molecular 
bands). In addition, they are very faint and therefore difficult to 
observe at high resolution and signal-to-noi se ratio (S /N) (e.g., 
IWoolf & Wallersteinll200HlBean et al.ll2006h . 

K-dwarf spectra are also useful for radial velocity determina- 
tions due to the large number of spectral lines present. Currently, 
radial velocity uncertainties of a fraction of 1 km s _1 are consid- 
ered reasonable for studies of kinematics and stellar abundances. 
Even with the highest S /N and highest resolving power spec- 
tra, these uncertainties cannot be reduced, in an absolute sense, 
below a level of a few hundreds of meters per second due to 
systematic errors. One of these is the relativistic gravitational 
redshift, which can be, however, reasonably well estimated from 
stella r mass and radius d eterminations, at least for nearby stars 
(e.g.. lDravins et alll!999l) . The other important source of uncer- 
tainty is due to granulation. 

Typically, granulation effects introduce a radial velocity un- 
certainty on the order of several hundreds of meters per second, 
an error that depends on the lines used, given that wavelength 
shifts depend on the properties of the spectral line, in partic- 
ular line strength. If both gravitational redshift and convective 
blueshift are properly taken into account, the uncertainty in the 
absolute radial velocity of the standard stars, some of which are 



K-dwarfs (e.g.. lStefanik et al.lll999l) . could be reduced dramati- 
cally to only a few tens of meters per second. 

Standard stars with very accurate absolute radial velocities 
can be crucial for the interpretation of data fro m massive kine- 
matic (and abunda nce) surveys suc h as RAVE (Steinmetz et al. 
2006), GAIA (e.g., Perrvman l l2005l), APOGEE dMajewski et al. 
120071) . and SPSS/SEGUE dRe Fiorentin et al.ll2007l) . In partic- 
ular, constant radial velocity objects can replace comparison 
lamps as wavelength calibrators, if needed (for example, this has 
been proposed for the GAIA mission, where the observations of 
a group of "well-behaved" s tars will be used to self-calibrate 
its very stable spectrograph; Kat z et al]|2004l) . The determina- 
tion of accurate radial velocities is also important for the de- 
tection of long term variability (e.g., unresolved wide binaries). 
Furthermore, the elimination of systematic errors from spectro- 
scopically determined radial velocities is necessary to assess the 
reliability of non-spectroscopic methods of radial velocity de- 
termination, such as those involving astrometric measurements, 
for example via secular parallax/pr oper motion variation s or the 
varying extent of moving clusters dDravins et al.ll 1999b . While 
these techniques were envisioned as early as the 1900's, their 
applicability has been limited by the accuracy of the astromet- 
ric da ta. Currently, for example with Hipparcos (Perrvma n et al.l 
1 1997b data, only the moving cluster method can in some 
cases achieve accur acies below 1 kms" 1 dLTndegren et al. 2000; 
iMadsen et al.l 12002) but the advent of future space astrometry 
missions (in particular GAIA), which will improve the accuracy 
of stellar parallaxes and proper motions by at least two orders of 
magnitude, guarantees that these techniques will become very 
powerful and the results widely used. 

In this series of papers, we perform a detailed study of gran- 
ulation in K-dwarfs and explore its impact on the observed spec- 
tra, in particular line shapes and wavelength shifts, as well as 
spectroscopic abundances and fundamental parameters. In this 
paper (Paper I), we look for the basic signatures of granula- 
tion on very high quality spec tra of a small sample of K-dwarfs. 
iPaper III (iRamfrez et al.l l2009al) deals with hydrodynamic sim- 
ulation s of K-dwarf envelopes, and in Paper III (Ramire z et al.l 
2009b) we explore the impact of 3D effects on the determination 
of stellar abundances and fundamental parameters. 

The data acquisition and processing are described in de- 
tail in th is paper. Th erefore, some sections are more relevant to 
Paper II or Pap er IlJ than to this paper. For example, the discus- 
sion about continuum normalization (Sect. |3.4b will be important 
for our stu dy of 3D effects on spectroscopic abundance determi- 
nations in Paper III while the details about the instrumental pro- 
file (Sect. [33} will be cr ucial for our comparison of predicted to 
observed line profiles in Pap er III 

2. Detection of granulation signatures in K-dwarf 
spectra 

At least 5 independent points across a spectr al line are re- 
quired to define its bisector (e.g jDravin s 1987b). Thus, for typ- 
ical strong lines in K-dwarf spectra with a full width of about 
10 kms -1 , we need at least 2 kms -1 resolution, or R = A/AA = 
150,000. In practice, however, the feasibility of measuring line 
bisectors also depends on the magnitude of the granulation ef- 
fects and the level of noise. In K-dwarfs, granulation effects 
are very small, on the order of tens to a few hundreds of me- 
ters per second. Numerical tests show that only at very high 
resolving power (R > 150,000) and high signal-to-noise ratio 
(S IN > 1 , 000) the line bisectors associated with granulation are 
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well defined. To reach such high S IN, however, very long ex- 
posure times are required for typical K-dwarfs (even the bright- 
est, nearest ones), which makes the observations impractical for 
mid-sized telescopes. At S /N - 300 and R ^ 200,000, the line 
bisectors are noisy, but the overall shapes of the bisectors due to 
granulation are still detectable. 

In addition to the limitations due to finite S IN, line blends, 
which are particularly important in K-dwarfs, can severely dis- 
tort the line bisectors. One way to overcome the effect of blends 
is to group bisectors of lines of the same species and similar 
strength, given that they form in roughly the same photospheric 
layers and are, therefore, expected to experience similar granula- 
tion temperature and velocity fluctuations. In this context, wide 
spectral coverage observations are ideal because they allow us to 
select a large number of clean lines from each spectrum. 

For the measurement of wavelength shifts, on the other 
hand, it is necessary to have, in addition to very well calibrated 
spectra, accurate laboratory measurements of rest wavelengths. 
Unfortunately, these are currently at the level of about 75 m s _I 
(Sect. 14.21 1. and our results are therefore subject to these uncer- 
tainties. The type of study presented here would therefore greatly 
benefit from efforts by laboratory physicists to improve these 
wavelength measurements. 



3. Observations 

3.1. Target selection 

The S pectroscopic Survey of Sta rs in the Solar Neighborhood 
(S 4 N, Ullende Prieto et al]|2004l) contains all G and early K- 
dwarfs (K2V and earlier) within 15 pc. Cooler stars are also 
included in the survey but their sampling is not complete. 
We used this catalog, and the stellar parameters derived there, 
to select our targets by constraining their effective tempera- 
tures to the range 4600 K < T e g < 5250 K and metallicities to 
-0.25 < [Fe/H] < +0.25, in addition to their observability from 
McDonald Observatory. Nine K-type dwarf stars were selected 
and are listed in Table Q] 

We also obtained spectra of Procyon and the day skylight, 
the latter as a proxy for the solar spectrum, to check the accu- 
racy of our data reduction procedures, given that these stars have 
very high quality spectra available from p revious studies (e.g., 
iGriffin & GnlM 19961: iKurucz et al.lll984k Note, however, that 
the skylight spectrum obtained in this manner may be slightly 
different than the real solar spectrum because the sky fills the 
spectrograph slit completely, therefore making the data reduc- 
tion different from that of a point source, and it is known that 
aerosol and Rayleigh-Brillouin sc attering affect s lightly the ob- 
served line shapes and strengths (Gray et al. 2000|). 

3.2. Observations and data reduction 

The data for our sample stars were acquired in 5 observing runs 
between February 2006 and May 200 7. The spectra w ere ob- 
tained with the 2dcoude spectrograph dTull et al.lll995l) on the 
2.7 m Harlan J. Smith Telescope at McDonald Observatory, us- 
ing the cs21 mode (grating e2, focus fl), which is expected to 
deliver a spectral resolution of R ^ 210,000 (see, however, 
Sect. 13.51 ). In the cs21 mode, a single exposure results in an 
echelle spectrum of about 20 orders, each of them covering ap- 
proximately 20 A. By rotating the spectrograph grating to 7 dif- 
ferent positions (hereafter referred to as "instrumental setups"), 
full coverage from 5580 to 7900 A was achieved. 



P revious observat i onal studies of granulation in K-dwarfs 
(e.g., lDravinsl[T987i lGraylll982l 120051) have used very high 
quality data but with an otherwise limited spectral coverage, 
sometimes involving the analysis of only one spectral feature. 
Although the quality of our observations is similar or slightly 
superior to that of these pioneering studies, our wavelength cov- 
erage surpasses all of them. 

A typical single exposure of 20 minutes for a V - 6 star re- 
sulted in S IN - 100 under good observing conditions. To reach 
our required signal-to-noise ratio (S IN > 300), several single ex- 
posures of the same setup were co-added before merging them 
into a single spectrum. Sects. [3~3l and l3.4l explain the details of 
the co-adding and merging procedures. 

Data reduction (bad pixel removal, overscan correction, flat- 
fielding, and scattered light correction) was performed using the 
standard IRAB3 echelle package. Comparison lamp (ThAr) 
spectra were obtained approximately every hour to ensure an 
accurate wavelength calibration (5 x 10~ 4 A of RMS scatter on 
average). Instead of linearizing the wavelength solutions (i.e., 
resampling each order to a wavelength scale of constant step), 
the actual pixel-to-wavelength relations were preserved to mini- 
mize interpolation errors in the observed counts. Fig.[T]shows an 
example of a reduced spectrum. 

In addition to the 2dcoude data, we obtained spectra o f our 
samp l e stars using the H igh Resolution Spectrograph (HRS; Tull 
1998; iHill etaT]|2006l) on the 9.2 m Hobby-Eberly Telescope 
(HET) for wavelength calibration purposes (details are given in 
Sect. 13.4b . This instrument is capable of delivering R ^ 120, 000 
spectra with complete wavelength coverage from 6000 to 7900 A 
(red chip) and from 4100 to 5900 A (blue chip) in one single ex- 
posure. The process of data reduction was similar to that applied 
to the 2dcoude data. For these spectra, the accuracy of the wave- 
length solutions is about 10~ 3 A for the red chip and 5 x 10~ 4 A 
for the blue chip. Exposure times were set to reach S IN - 750 at 
the center of the red chip. The procedures described in the next 
section were also applied to these data. 

3.3. Shifting and co-adding: increasing the SIN ratios 

As explained above, several individual exposures of each setup 
were required to reach high signal-to-noise ratios. Co-adding 
them is not trivial because even after correcting for the Earth's 
motion, wavelength shifts due to the intrinsic radial velocity of 
the star (e.g., due to binarity, variability, and presence of planets) 
or instrumental and environmental effects (the 2dcoude spectro- 
graph is not inside a vacuum chamber) are present between expo- 
sures. These shifts are on the order of tens to hundreds of meters 
per second and, therefore, of the same magnitude as the effects 
we are trying to resolve. Consequently, these shifts must be taken 
into account before co-adding in order to keep the spectra useful 
for granulation studies. 

We determined radial velocity shifts between pairs of spectra 
of the same object by cross-correlation in the Fourier space, af- 
ter rebinning the spectra to have a constant step in log A. Orders 
severely affected by telluric lines (more than 50% of the A range) 
were discarded. Our cro ss-correlation calcu lations were made 
using the IDL xc code bv lAllende Prietol (120071) and not the pop- 
ular IRAF fxcor task because the latter requires the spectra to 
be linearized, a procedure that we avoided (see Sect. 13.21 ). 

1 IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with the 
National Science Foundation - http : / /iraf . noao . edu 
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4889 
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6.4 
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HD 160346, GJ688 
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HD 219134, HR 8832, GJ 892 
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K3 


6.5 


4641 
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Table 1. Sample observed in this work. Basic data (V, SpT) were obtained from SIMBAD. Distances and absolute magnitudes are 
based on Hipparcos parallaxes. The atmospheric parameters r e ff and [Fe/H] are those determined bv lAllende Prieto et alj (|2004) 
and have typical uncertainties of about 100 K and 0.1 dex, respectively. 
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Fig. 1. Sample spectrum from our 2dcoude-cs21 observations for HIP 96100. Only one instrumental setup is shown. 



Determining the uncertainty of our cross-correlation proce- 
dure is a difficult task. Monte Carlo simulations using both syn- 
thetic spectra and an observed high S IN spectrum predicted er- 
rors of only about 4 m s" 1 for R = 200, 000 and S /N = 100. This 
procedure leads to optimistic error estimates because it does not 
take into account other effects such as pixel shifts that are non- 
linear in velocity and order-to-order scatter due to instrumental 
distortions, the fact that spectral lines are not symmetric or at 
their rest wavelengths due to granulation effects, incompleteness 
or uncertainties in the atomic data used for the calculation of the 
synthetic spectra, telluric lines, etc. Instead, we estimated the er- 
ror from the order-to-order scatter only, for which we found no 
strong correlation with wavelength. Thus we find that our rela- 
tive radial velocities have a mean error of about 14 m s~'. 



The standard approach to the determination of radial ve- 
locity shifts between spectra is to cross-correlate them with a 
reference spectrum of known radial velocity. In our case, we 
do not have access to a reference spectrum of higher quality 
than the data (e.g., high resolution synthetic spectra including 
line asymmetries over the whole spectral range of our observa- 
tions) and therefore we need to adopt a different approach. In 
this work, the number of spectra co-added for a single faint ob- 
ject in a given setup was typically between 5 and 10. Instead of 
cross-correlating each of them with a single reference spectrum, 
we cross-correlated all of them with each other. We then used 
all the measurements to refine the shifts of all spectra with re- 
spect to the first one of the list employing "self-improvement" 
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dAllende P rieto 2007). Using this technique, the mean error in 
the relative radial velocities was reduced by about 15%. 

After correcting for the radial velocity shifts, the spectra 
were coadded using an algorithm that minimizes interpolation 
errors in the observed fluxes. When coadding, each original pixel 
has to be interpolated to some common wavelength scale. One 
may be tempted to use the wavelength solution of one of the 
spectra and interpolate the rest to that dispersion at the risk of 
introducing large interpolation errors. Instead, we determined, 
for each set of spectra to be coadded, a wavelength solution such 
that the flux interpolation distances are minimized. 

The pixel-to-wavelength relation of the first spectrum was 
used as a starting reference. The mean difference between adja- 
cent pixels, 6A, was then determined. For each pixel in the ref- 
erence spectrum, we defined a box extending to +6 A/2 around 
its wavelength. Next, we looked at the rest of the spectra and 
found all pixels with wavelengths within that box. The average 
wavelength of all those pixels was then adopted for that particu- 
lar position. In this way, the mean interpolation distance was re- 
duced by a factor of 2, compared to the case in which the spectra 
are interpolated to the dispersion of the first (or any other ref- 
erence) spectrum. The interpolation distance for the linear case 
(i.e., rebinning the spectra to a common constant step disper- 
sion) is substantially larger (a typical factor of 5 to 10) but, as 
explained above, we avoided this type of interpolation. 

3.4. Normalization and "merging" 

As explained in Sect. 13. 1| all of our sample s tars are included in 
the S 4 N survey of Alle nde Prieto et all d2004l) . which provides 
atlases of spectra covering the wavelength range from 362 to 
921 nm with a resolving power R 60,000. The continuum 
normalization of these data is superb; not only were the blaze 
shapes removed by fitting high order polynomials to the up- 
per envelopes of the observed fluxes in each order, but also the 
smooth variation of the blaze shapes in the direction perpendic- 
ular to the dispersion was taken into account, thus removing the 
instrumental response with a very careful tw o-dimensional mod- 
eling, as described in lBarklem et al.l d2002). This normalization 
was essentially inherited by our spectra using the procedure de- 
scribed below. 

Each order of our coadded cs21 spectra (one instrumental 
setup at a time) was divided by its corresponding cs23 piece from 
the S 4 N data set, after lowering the resolution of the cs21 data to 
match that of the cs23 data, correcting for radial velocity shifts, 
and rebinning to a common wavelength scale. In principle, the 
result should be a smooth function corresponding to the shape of 
the continuum (the upper "envelope") of the cs21 observations. 
However, due to the finite S /N values and temporal variations of 
the strengths of some lines (in particular the telluric ones), this 
envelope had to be smoothed out using a median filter of width 
equal to 100 pixels (about 2 A). The cs21 data were then divided 
by this envelope in each order, for each instrumental setup. 

The coadded and normalized cs21 data (7 setups for each 
object) were finally merged using an HET spectrum of the same 
object as a radial velocity template. Note that the spectra for each 
setup can have a different radial velocity, since we only removed 
relative, and not absolute, Doppler shifts (Sect. [331 ). Therefore, 
each order of a cs21 setup was cross-correlated with its corre- 
sponding piece in the HET-HRS spectrum, after lowering the 
resolution of the former to match the R ~ 120, 000 of the HRS, to 
determine the radial velocity needed to align the cs21 spectrum 
with the HET-HRS spectrum, which has a unique wavelength 



solution. The robust (iterated and weighted) average shift of all 
orders for each setup was adopted to correct the cs21 data, thus 
removing inaccurate shifts from orders affected by telluric lines 
or instrumental defects. In the regions of overlap, a weighted 
mean for the fluxes was calculated after adopting the common 
robust rebinning explained in Sect. 13. 31 Fig.|2]illustrates this step 
of the data processing. 

On average, the order-to-order scatter of the radial velocities 
used for the merging of setups was 70 m s _1 . This error is larger 
than the 12 m s error quoted in Sect. |3.3| for the order-to-order 
scatter of the cs21 data alone, in part due to the lower resolving 
power of HRS, but also probably because of distortions in the 
wavelength-to-pixel mapping of HRS. Note that the regions of 
overlap of the different cs21 setups are small compared to the 
wavelength coverage of each order. Only within the regions of 
overlap, spectral line shapes are affected by this large merging 
error of 70 m s _1 . The majority of spectral lines (approximately 
80% of them), however, are not in the regions of overlap and 
the wavelength mapping within a spectral line is still accurate 
at the 12 m s _1 level. On the other hand, in the final merged 
spectrum, the absolute wavelength scale of all setups and orders 
combined is uncertain to about 70 m s _1 . This 70 m s _1 uncer- 
tainty will have an effect on our measurements of wavelength 
shifts (Sect. l4~2l but not on the line profiles or their bisectors. 

3.5. Resolving power and instrumental profile 

To calculate the resolving power of 2dcoude-cs21, we used the 
ThAr exposures. We measured the full width at half maximum 
(FWHM) of indi vidual Th ( fai ) and A r (/at) lines using the line 
identifications of Murph y et al.l(l2007l) . Note that the Th lines are 
sharper than the Ar lines due to their different atomic weights, 
which makes thermal broadening more important for the lat- 
ter. The thermal broadening (cr) is inversely proportional to the 
square root of the atomic weight A; thus 



OTh 
CAr 



Am 



1/2 



= 0.415. 



If the instrumental FWHM is /, then 

/ + cr Th = / Th 



/Ar 



(1) 



(2) 
(3) 



from which we obtain an expression for the instrumental FWHM 
that is independent of temperature: 



f 2 ■ 



0.17/* 



0.83 



(4) 



The resolving power is then A/f. 

The measured FWHM values were found to vary signifi- 
cantly with time, even within a given run of a few days. From 
one run to another (months apart), the change in the average 
FWHM values was significant, as shown in the upper panel of 
Fig. [3] Temperature broadening effects alone cannot explain the 
diversity of values seen in Fig.[3]given that less than 5% of the 
broadening of Th lines is thermal (see lower panel of Fig. [3). 
Considering changes of 10% in the temperature of the ThAr 
lamp (which was actually well controlled), the total broadening 
would be affected by only 0.05 x V0.1 = 1.6% or less, yet in Fig. [3] 
(upper panel) we see variations of up to 20%. The minimum and 
maximum values of the resolving power of our observations are 
about 160,000 and 210,000, respectively. 
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Fig. 2. Merging of spectra from different instrumental setups (Pl-7). Only one order per setup is shown here; note that each setup 
has about 20 orders in different wavelength regions (see Fig. |T}. The fluxes have been shifted arbitrarily in the vertical direction for 
clarity. The cs21 pieces (Pl-7) have been aligned with the HET+HRS spectrum by means of cross-correlation and are ready to be 
merged. 
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Fig. 3. Upper panel: Average FWHM of the instrumental profile 
for each of our observing runs. The 2dcoude spectrograph (cs21 
mode) is expected to deliver / ^ 1.4 kms -1 at all times; the 
differences shown here reveal that the spectrograph has suffered 
from focus degradation. Lower panel: Average ratio of instru- 
mental to Th line FWHM values for each of our observing runs. 
The constancy of this ratio suggests that the differences seen in 
the upper panel of this figure are not related to temperature ef- 
fects on the ThAr lamp. 



Note that in Eqs.|2]and[3] as well as in the discussion given 
above, it was assumed that the instrumental profile of the spec- 
trograph was Gaussian. Although the actual ThAr line profiles 



are slightly asymmetric, Gaussian fits made reasonably good ap- 
proximations, as shown in Fig. 2] where the average Th line pro- 
files for each observing run are shown. They were obtained from 
the individual Th line profiles, normalized to 1 at their maxi- 
mum before averaging in bins of 0.5 km s _1 , which corresponds 
to about one-third of the width of the resolution element. The de- 
parture from Gaussianity of the real instrumental profile shown 
in Fig. [4] is very small. The most prominent non-Gaussian fea- 
ture is observed on the blue wing of the instrumental profile (near 
-2 km s~') and it amounts to less than 10 %. In fact, when the 
instrumental profile was sharpest (in February 2006), the maxi- 
mum difference between the real profile and the Gaussian profile 
was only about 5%. A similar bump on the red wing is observed 
only for the May 2007 observations, which correspond to the 
most degraded (wider) instrumental profile of our observations. 
This information will be crucial for the comparison of observed 
line profil es with th ose predicted by hydrodynamic model atmo- 
spheres in Pap er 11 



4. Measurement of the granulation signatures 

The spectra of K-dwarfs, in particular those with near solar and 
super-solar metallicities, are very rich in absorption lines. For 
spectroscopic studies of granulation, iron is an ideal element; it 
has a large atomic mass, which implies low thermal broaden- 
ing and therefore a high sensitivity to non-thermal fields, one 
isotope is much more abundant than the others, thus reducing 
the impact of asymmetries due to fine structure, it is abundant, 
and its atomic structure is very complex, allowing a very large 
number of transitions. Furthermore, contrary to most other heavy 
elements, reliable laboratory data are available for iron. In this 
paper, we deal with the detection of the granulation signatures 
using Fe i lines exclusively. Although the granulation effects are 
expected to be stronger in Fe n lines owing to their larger for- 
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Fig. 4. The average instrumental profile for each of our observ- 
ing runs is shown with diamonds connected by dotted lines. The 
solid lines correspond to Gaussian profiles of FWHM identical 
to the average FWHM of the instrumental profiles, which are 
given in the upper right corner of each panel in km s _1 . 

mation depth, few Fe n features are available in our spectra and 
almost all of them are severely affected by blends. 

4.1. Line bisectors 

We measure d the bisectors of a large number of spectral lines 
listed in the iNave et all (1 19941) multiplet table for Fe i. In our 



available spectral range (5580 < A < 7800 A). INave et al.l list 
about 900 lines, most of which are weak or blended. For each 
spectrum, we examined the neighboring ± 1 A of all ~ 900 Fe i 
lines mentioned above and attempted to measure their line bi- 
sectors. Very weak lines with minimum normalized flux values 
above 0.98 were discarded because their bisectors are very sen- 
sitive to blends and noise. Spectral lines in regions where the 
local S IN was below 200 (as estimated from the total number of 
counts in the local continuum) were also discarded, except for 
HIP 37349, which has S/N 150 everywhere. 

The line bisectors were measured by determining iso-flux 
wavelength (or velocity) points on each wing of the line, re- 
sampling the original profiles using cubic spline interpolations. 
The error bars were obtained by propagating the observational 
erro rs of our sp ectra, as quantified by the local continuum S /N 
(see lGravTl 983 for details). Hereafter, the line bisectors are mea- 
sured with respect to the line core wavelength. The latter mea- 
surements are discussed separately in Sect. 14.21 

Many of the Fe i lines given in the INave et all catalog are 
blended in the observed stellar spectra. In many cases the blends 
will be clearly distinguishable by visual inspection and a line 
selection could be made in this manner. However, even weak 
line blends can significantly distort the shape of the line bi- 
sectors, and, given that the line asymmetries are very small in 
K-dwarfs, these blends should be avoided using a different ap- 
proach. Instead of selecting subjectively spectral lines based on 
visual inspection of the line itself or its bisector, we adopted an 
automated scheme in which the smoothness of the line-bisector 
is used as evidence that the line is not blended. For each bisec- 
tor measured, we computed the slope d{A\)/df n , with Av be- 
ing the abscissa values of the line bisector and /„ their corre- 
sponding normalized fluxes. Bisectors for which d(Av)/df n > 
4, 000 m s _1 per normalized flux unit were discardedQOf course, 
lines blended with features that share the same core wavelength 
will survive this criterion but their number will be small com- 
pared to blends on the line wings and, furthermore, even sub- 
jective criteria will have difficulty disentan gling them from true 
clean lines. About 40% of the lines listed by Nav e et all survived 
our selection criterion. 

Even after this selection of smooth bisectors, the level of 
noise was relatively high, as expected for the typical signal-to- 
noise ratios of our observations. Since the shape of the line bi- 
sector depends primarily on the line strength, the noise can be re- 
duced by averaging bisectors of spectral lines of similar strength. 
This is shown in Fig. [5] for the case of the star HIP 96100. 
Although a wide range of bisector shapes is apparent in Fig. [5] 
we applied robust statistics to obtain the mean bisectors, itera- 
tively discarding points that disagree from the mean values by 
more than 2.5<x, which makes our final results reliable. Similar 
mean Fe i line bisectors were computed for all the K-dwarfs in 
our sample and are shown in Fig. [6] 

It is reasonable, although not strictly correct, to assume that 
the line bisectors of lines of similar strength are identical. In 
that case, the standard deviation (<x) is simply a measure of the 
noise and the actual measurement error of the mean line bisec- 
tors should be given as the standard error (<x/ y/nj) of the n; line 
bisectors used in each group. In Fig. the error bars correspond 
to the standard error. 



2 On average, we find that bisectors of strong lines have d(Av)/df„ ^ 
150 m s -1 per normalized flux unit and therefore this criterion is a con- 
servative but robust one. 



8 



Ramirez et al,: Granulation in K-Dwarfs. I. 




-400 -200 



200 400 



Velocity (m s ) 




-400 -200 



Velocity (m s ) 




HIP 96100 



200 400 



-400 -200 



200 400 



Velocity (m s ] ) 



Fig. 5. Mean bisectors for HIP 96100 (black solid lines with error bars) obtained from individual line bisectors (gray lines). The 
error bars correspond to a Icr scatter. 



The mean line bisectors for HIP 96100 have the largest 
spar0 among our sample stars. There are three reasons that can 
explain this. First, the star has the hottest effective tempera- 
ture in the sample, second, its metallicity is lower than solar 
([Fe/H] = -0.22), and, finally, its projected rotational veloc- 
ity is the lowest among our sample stars (V sin / ^ 0.8 km s _1 , 
see Appendix [At. The granulation effects, and in fact the abso- 
lute magnitude of the inhomogeneities, are predicted to increase 
with effective temperature and also with decreasing metallicity, 
although probably the la tter effect is significantly smaller (e.g., 
lAllende Prieto et alj 19991) . while low projected rotational veloc- 
ities ( V sin i < 1 km s _1 ) do not alter t he shape of the line bisec- 
tors significantly, as will be shown in Paper II. The dependence 
on e ffective tem perature has in fact been observed in real stars 
(e.g.. lGravl!99a> . 

In this context, it is interesting to see the case of the next 
few hotter stars, HIP 26779, HIP 16537, HIP 88601, HIP 64797, 
and HIP 37349, which do not show bisectors with large spans. 
In fact, in some cases they are smaller than those of the coolest 
star in our sample (HIP 2331 1). The reason for this is probably 
their large V sin i values (approximately between 2 and 3 km s _1 , 
see Appendix lAb . which reduces the span of the line bisec- 
tors (Paper II). We note here that the data for HIP 64797 and 
HIP 37349 are noisier than the average. Most of the observa- 
tions for HIP 64797 were made under poor weather conditions 
while the spectrum of HIP 37349 has only S/N ~ 150. 

In addition to the stellar rotation effect in the 5 stars men- 
tioned above, it is also interesting to note that they show high 
levels of chromospheric activity (with the possible exception 
of HIP 88601), as evidenced by their Ca n H and K line pro- 
files (Appendix iBt. Chromospheric effects on spectral line for- 
mation have not yet been studied in enough detail to deter- 
mine quantitatively their impact on our line bisector measure- 
ments. Nevertheless, the fact that the signatures of granulation 
seem to be slightly different in the most active stars of our sam- 
ple suggests a correlation. An active chromosphere may have 
an impact on the line shapes and intensities in several ways. 
Emission may fill the cores of some lines, as it is obvious in 
the H and K lines, but it may be also present to a lesser ex- 
tent in weaker lines. Zeeman splitting could be m easurable due 
to the presen ce of strong magnetic fields (e.g., iGravl Il984bt 
Borrero 2008), while magnetic pressure may have an impact on 
the photospheric structure. We can conclude that the line bisec- 
tors are significantly affected by the presence of strong mag- 
netic fields and/or stellar activity, although the magnitude and 



1.0 

0.8 
0.6 
0.4 
0.2 
1.0 
0.8 
0.6 
0.4 
0.2 




skylight 




HIP 37279 



-400 -300 -200 -100 
Velocity (m s -1 ) 



3 We define the "span" as the difference in velocity between the red- 
dest and bluest points of the line bisector. 



Fig. 7. As in Fig. |6]for our observations of the Sun and Procyon. 



impact of these effects are largely unknown for stars other than 
the S un, where this effect has been documented in some detail 
(e.g., Balasubramaniam 2002i [Schlichenmaier et al.l l2004). 

HIP 86400 is our reference star because it has parameters 
very close to those adopted in t he calc ulation of the 3D model 
atmosphere that we will study in lPaper iH While the mean bisec- 
tor for the weak lines shows no significant asymmetry (a span of 
only about 15 m s ), the strongest lines show bisectors that 
span about 75 m s on average. 

The coolest stars in our sample, HIP 1 14622 and HIP 2331 1, 
both have relatively low V sin i (< 2 km s -1 , AppendixlAl. Their 
spectra are richer in spectral lines compared to those of our other 
sample stars and this is perhaps the reason why their mean bisec- 
tors show large deviations from the expected behavior near the 
continuum. In fact, the bisectors approach zero velocity there, 
which points to the influence of numerous blends randomly dis- 
tributed on both sides of the spectral lines. Nonetheless, far from 
the continuum, their line bisectors are similar to those of the ref- 
erence star HIP 86400. Therefore, it is likely that the granulation 
contrast does not change significantly in dwarf stars of effective 
temperature between 4600 and 4800 K. 

To check that our measurements of line asymmetries are not 
affected by observational or data processing artifacts, we ob- 
tained mean line bisectors in our skylight and Procyon spectra 
as well (Fig. [7]). The mean bisectors of the skylight are much 



Ramirez et al.: Granulation in K-Dwarfs. I. 



1.0 



0.6 ■ 



en 0.4 - 



0.2 ■ 




HIP 96100 



H+ 



+++ 



ft 





HIP 26779 
I I + ' 



I I 




HIP 16537 



1.0 ■ 



0.6 ■ 




01 0.4 ■ 



0.2 - 



HIP 88601 
I I " 



I l^-J 







HIP 64797 

I I ) I' 



HIP 37349 




1.0 ■ 



0.8 ■ 



0.6 ■ 



cn 0.4 



0.2 ■ 







HIP 86400 



HIP 1 14622 



HIP 2331 1 




-400 -300 -200 -100 
Velocity (m s -1 ) 



-400 -300 -200 -100 
Velocity (m s _1 ) 



-400 -300 -200 -100 
Velocity (m s -1 ) 



Fig. 6. Mean line bisectors for our sample of K-dwarfs, measured with respect to their core wavelengths. Shifts of -200 and 
-400 m s _] have been applied to the weaker lines for clarity. The error bars correspond to the standard error of the line-to-line 
scatter, i.e., crj -\fn~i, where cr is the standard deviation and n/ the number of lines. Bisector points for which <x > 200 m s _1 have 
been discarded (these often appear near the continuum). 



larger than those of our reference star. Note, however, that our 
skylight observations are significantly different than stellar ones. 
Solar line bisectors have been determ ined using the actual so- 
lar spectrum ( e.g., Kurucz et al. 19841) by several authors (e.g., 
Dravins et al.1 1 198 it lAsplund et alJ 120001 lAllende Prieto et all 



2002), who show that their amplitudes are about twice as large 



as those shown in Fig. [7] although this is in part due to the fact 
that their observations h ave higher spectral resolution (for ex- 
ample, the lKurucz et al.l atlas has a resolving power of 500,000). 
The mean bisectors measured in our Procyon spectrum, however, 
are in excellen t agreement with those measured elsew here (e.g., 
lDravinsll987atlGravll989tlAllende Prieto et al.l2002l) . Note that 
the metallic lines in the spectrum of this hotter star are, in gen- 
eral, weaker than those observed in K-dwarfs and, in our spec- 
tral range, very few strong lines are available. This is the reason 
why Fig. [7] does not show the strongest mean line bisector for 
Procyon. 



4.2. Wavelength shifts 

The second fundamental test of granulation signatures in stel- 
lar spectra involves measuring core wavelength shifts. A large 
convective blueshift is expected to occur for the weakest lines, 
which are formed in deep layers and see strong granulation con- 
trasts, while the cores of strong lines are much less affected given 
that they form in the highest photospheric layers, where the cor- 
relation between intensity and velocity fields is weak. We should 
therefore look at the relation between line shift and line strength 
(quantified, for example, by the line equivalent width) to deter- 
mine the magnitude of the granulation effects and compare them 
to the model predictions. 

Unfortunately, there are large uncertainties in the determina- 
tion of the line shifts from our observed spectra, some of which 
are independent of the quality of our data. For example, very ac- 
curate laboratory wavelengths for l arge numbers of sp ectral fea- 
tures are necessary. The catalog bv lNave et al.l (fl994) is perhaps 
one of the most accurate and complete line lists currently avail- 
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Fig. 8. Core wavelength shifts measured in the spectra of our sample stars as a function of line equivalent width (small diamonds). 
The filled circles with error bars represent average values in bins of equivalent width. The zero point of this relation has been set 
using the strongest lines, those that define the plateau. 



able for this purpose. Nevertheless, the wavelengths measured 
by them, in the best cases, have uncertainties of about 75 m s _1 . 

We measured the c ore wavelengths of as many Fe i lines as 
possible included in the lNave et alj d 19941) catalog and restricted 
our measurements to those lines that have wavelengths measured 
with uncertainties lower than about 75 m s . The core wave- 
lengths were obtained by fitting fourth order polynomials to the 
1 1 points closest to the flux minimum. A more detailed descrip- 
tion of the method we adop ted to determine core wa velengths 
and their errors is given bv lAllende Prieto eta il (120021) . 

Since the wavelength shifts i n the solar spectrum have been 
measu red with accuracy (e.g., lAllende Prieto & Garcia Lopezl 
1998a) and the convective blueshift vs. equivalent width rela- 
tion is well defined there, we used our results for the skylight 
spectrum to perform a line selection. Only those lines for which 
the convective blueshifts were between 2cr of the mean line shift 
vs. equivalent width trend were used. Statistically, this selection 
reduces the impact of line blends and misidentifications. 

Our measured line shifts could be systematically displaced 
with respect to their absolute values due to the gravitational red- 
shift and the uncertainty in the absolute stellar radial velocity. 
Thus, we adopted an arbitrary zero point using the core wave- 



length shifts of the strongest lines. Basically, the mean line shift 
of the strongest lines was set to zero (see below for details). 

The line shift vs. equivalent width relations for our sample 
stars are shown in Figs.[8]and[9] In each panel, we show the indi- 
vidual line shifts as well as averaged values in bins of equivalent 
width (the filled circles with error bars) to allow a better visual- 
ization given the large scatter observed. The bins have a width 
of 10 mA below EW = 100 mA and 35 mA above it. When the 
number of data points per bin is lower than 3, the bin-averaged 
values are not shown. 

The basic signature of granulation is evident in all cases. The 
strongest lines show the smallest blueshifts while weaker lines 
are significantly blueshifted. Furthermore, the slope in the line 
shift vs. equivalent width relation is larger for the skylight com- 
pared to the K-dwarf observations and even larger for Procyon. 
Although the scatter for the latter is large and there is not enough 
line strength coverage to draw in reasonably well the actual form 
and zero point of the trend (Fig.|9]l, this has been established by 
previous studies (e.g. . lAllende Prieto et a l. 2002). 

In several cases (e.g., HIP 96100, HIP 26779, HIP 86400, 
HIP 2331 1) it is clear that the trend is nearly linear for the weak 
lines but it reaches a saturation point at about 100 mA, from 
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Fig. 9. As in Fig. [8]for our observations of the Sun and Procyon. 



where it becomes a plateau. The detection of this plateau is of 
utmost importance, as it will be argued in the next section. As we 
mentioned above, the line shifts were displaced arbitrarily so that 
the strongest lines have zero shifts. We define the set of strong 
lines as those belonging to the plateau, where it is detectable. 
In other cases we assume that lines with EW > 150 mA are 
"strong" in this context. 

For the reference star, HIP 86400, the convective blueshifts 
are about -150 m s _1 for lines of EW ^ 20 mA. The blueshifts 
decrease (i.e., the line shift becomes less negative) as we ap- 
proach EW ^100 mA. Beyond this EW value, the shift is inde- 
pendent of EW. 

The plateau is not clearly detected in HIP 16537, HIP 88601, 
HIP 64797, and HIP 37349. The latter is likely due to the low 
quality of the data compared to the rest of the sample stars. For 
the other stars in this list, this non-detection suggests that effects 
other than granulation shape the line profiles. Interestingly, ex- 
cept for HIP 26779, these are the same stars whose line bisectors 
have properties that differ from the theoretical expectation. As 
we suggested in Sect. 14.11 the fact that these objects are also the 
most active among our sample stars may point to a correlation. 

In setting the zero point of the convective blueshifts, we have 
assumed that the wavelengths of the cores of very strong Fe i 
lines in K-dwarfs are not shifted. This is a reasonable assump- 
tion considering that these lines form in very high photospheric 
layers where the correlation between intensity and velocity fields 
is weak. Moreover, the assumption is justified by the fact that, in 
the solar case, where core wavelengths can be measured on an 
absolute scale, t he plateau of convective blueshifts is f ound to be 
at zero velocity (Allende Prieto & Garcia Lopez 1998b). 

Although there is large scatter for small EW values, it ap- 
pears that the line shifts become less dependent on EW below 
EW ^ 40 mA, suggesting another plateau (this is clearly visible 
in the case of HIP 96100). Interestingly, this result is somewhat 
consiste nt with our 3D model predictions, as it will be shown in 
iPaper 11 



5. Impact on radial velocity determinations 

Although ultra high precision (close to the meter per sec- 
ond accuracy) can be achieved for relative measurements (e.g., 
iMarcv & Butlenll998tlOue~ioz et al.ll200lHEndl et al.ll20Q5h . ab- 
solute radial velocities are subject to a series of systematic er- 
rors that are very difficult to remove dLTndegren et al. 1999; 
iLindegren & Dravinsl2003l) . The first of these is the gravitational 
redshift, which is on the order of 0.6 km s _I in FGK stars. A 
secondary relativistic correction for the transversal Doppler ef- 
fect is also requ ired for high velocity stars, where it amounts to 
about 50 m s _1 dLindegren et al.l ll999). Both of these effects in- 
troduce a constant shift in the radial velocity zero point of the 
star. Equally or more important, however, is the effect of granu- 
lation. 

If the cross-correlation method is used to determine the ra- 
dial velocity of a star, then only when a very similar object is 
used as a radial velocity standard will the method be success- 
ful. This is due to the fact that the spectra of stars with differ- 
ent fundamental parameters show different line strengths, in ad- 
dition to suffering from different granulation effects because of 
their different intensity and velocity field contrasts (Nordlund & 
Dravins 1990), although the impact of the latter may be in gen- 
eral smaller. On the other hand, measuring the core wavelengths 
of a multitude of spectral lines and applying the Doppler formula 
to determine the stellar radial velocity introduces the error of the 
unknown convective blueshift to the measurements. However, 
as explained in Sect. 14.21 the cores of very strong Fe i lines in 
K-dwarfs have zero convective blueshift and, therefore, one can 
avoid the effects of granulation in radial velocity measurements 
by only using spectral lines whose convective line-shifts fall in 
the "plateau" region of the line-shift vs. equivalent width rela- 
tion, provided the plateau is reliably detected, as it is the case 
of some of our K-dwarfs: HIP 96100, HIP 86400, HIP 23311, 
HIP 88601, and HIP 26779 (cf. Fig.©. 

The average values of the radial velocities we determined by 
cross-correlation of our cs21 spectra with each other and with 
their corresponding HET spectrum (see Sect. [3] for details) are 
given in Table [2] While our relative velocities are very accurate 
(the relative velocities of cs21 pairs have a mean uncertainty of 
12 m s _1 while that of a cs21-HET pair is about 70 m s~'), the 
HET reference spectrum has been only approximately corrected 
for the stellar radial velocity by using the Doppler shifts mea- 
sured for a large number of strong spectral lines. The granula- 
tion effects are not yet removed from these data. Therefore, the 
radial velocities given in Table [2] are only first order approxima- 
tions, useful to check for large velocity variations, but should not 
be used as the radial velocities of the stars. 

Table|2]shows that almost all of our sample stars have nearly 
constant radial velocity, which makes them useful as radial ve- 
locity standards. HIP 88601, however, is a very long period vi- 
sual binary star, while HIP 26779 is a suspected long period ra- 
dial velocity variable (see details below). Only HIP 86400, our 
reference star, shows large, relatively short-term, radial velocity 
variations. Its orbital parameters are derived in AppendixICl 

We determined absolute radial velocities for the non-variable 
stars HIP 96100 and HIP 23311 using measurements of core 
wavelength shifts of their plateau lines. We also determined ab- 
solute radial velocities for HIP 26779, a suspected radial veloc- 
ity variable, and HIP 88601, which is the primary star of a visual 
binary system that has an orbital period of 88.3 years an d a ve- 
locity semi-amplitude of about 3.4 km s" 1 dHeintzll 1988b . which 
can be useful to improve the orbital solution of the systems or 
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Object v' r <r(v' r ) n 

km s" 1 



HIP 96100 


27.1 


0.1 


31 


HIP 86400 


23.3 


4.2 


72 


HIP 233 11 


21.9 


0.1 


41 


HIP 16537 


16.4 


0.3 


30 


HIP 64797 


7.9 


0.3 


106 


HIP 26779 


1.7 


0.2 


53 


HIP 88601 


-5.4 


0.2 


31 


HIP 1 14622 


-18.0 


0.1 


30 


HIP 37349 


-18.0 


0.2 


35 



Table 2. Approximate values of the radial velocity obtained 
for our sample stars using our data. The last column corre- 
sponds to the total number of observations. HIP 86400 is a con- 
firmed single-lined spectroscopic binary; its orbital parameters 
are given in Table IC.ll These values are not intended to repre- 
sent the most accurate radial velocity determinations (those are 
given in Table 0; they only serve our purpose of looking for 
obvious short-term radial velocity variations. 



look for second order variations in the future. The procedure is 
described in detail below. 

As explained in Sect. [3] one HET spectrum was used as a 
radial velocity template on which the cs21 pieces of the spec- 
trum were put together. This HET spectrum was first shifted to 
an approximate laboratory wavelength frame by using a veloc- 
ity obtained by comparing the observed core wavelengths of a 
large set of atomic lines with those measured in the laboratory. 
At that point, we were not interested in the absolute radial ve- 
locity of the star but only required to have the spectrum within 
about 1 km s _1 of its actual radial velocity for the purposes of un- 
equivocally identifying strong spectral features. We will denote 
this first approximation to the radial velocity of the star with 
Av™ T (see Table |3). The barycentric correction for the Earth- 
Moon system orbit around the Sun (Avj! ary ) was then determined 
using the method outlined by McCarthy] (119951) . The maximum 
error for this correction is expected to be about 5 ms 4 (M. Endl; 
private communication). 

In Fig. [8] the plateau lines have been put to zero velocity 
(on their average) by applying the AVr' at corrections given in 
Table [3] Since, in principle, the plateau lines have zero con- 
vective blueshifts, this correction is precisely what is needed to 
convert Av™ T + Aw ^ into the absolute radial velocity of the 
star, but still uncorrected for relativistic effects (yf). Note that 

the standard deviation <x from the mean estimate of Av^' 3 ' is be- 
tween about 30 and 60 m s -1 . Assuming that the line shifts of 
very strong lines in K-dwarfs are not affected by any other sys- 
tematic error, the standard error cr/ yfiT[ would be an appropriate 
estimate of the uncertainty of this correction (ni being the num- 
ber of strong lines used to determine the average value). Thus, 
after including the uncertainty of the barycentric correction, the 
error in our estimate of \f is only about 10 ms -1 . 

To finally determine the absolute radial velocity of these 
stars, the relativistic gravitational redshift must be removed. This 
correction is given by Av^ 1 = -cAA/A = [-636(M//?) + 3] m s -1 , 
if M and R are given in solar units. The +3 ter m corrects for 
the ef fects of the Earth's gravitational field (e.g., iDravins et alj 
119991) . We estimated the masses and radii of the four stars 
listed i n Table[3]from the theoretical iso chrones b y Bertelli "et al.l 
(fl994T). using the methods described in lReddv et al.l d2003l) and 
Alle nde Prieto et aT] d2004l) . Using the mass-to-radius ratios in- 
ferred from these calculations, we determined the Av^ el values, 



HIP 


JD-2454000 


Av het 


Av r baly 


AvP' at 


(std. error) 


96100 


231.936693 


25.45 


1.56 


-0.07 : 


± 0.04(0.01) 


23311 


406.891895 


7.58 


14.35 


-0.01 : 


± 0.03(0.01) 


88601 


205.922914 


-29.54 


24.81 


-0.03 : 


±0.06(0.01) 


26779 


350.969887 


-23.79 


25.19 


-0.15: 


±0.06(0.01) 



HIP 








■l 


Vr 




96100 


26.93 ± 


0.01 


-0.68 ± 


0.06 


26.25 ± 


0.06 


23311 


21.92 ± 


0.01 


-0.71 ± 


0.05 


21.21 ± 


0.05 


88601 


-4.76 ± 


0.01 


-0.65 ± 


0.06 


-5.41 ± 


0.06 


26779 


1.25 ± 


0.01 


-0.71 ± 


0.06 


0.54 ± 


0.06 



Table 3. Radial velocities of four of our sample stars. The several 
corrections are explained in detail in the text. The last column in 
the bottom table gives our best estimate of their absolute radial 
velocities. All velocities are given in km s _1 . JD is the barycen- 
tric Julian date. 



which allowed us to determine the absolute radial velocities, v f , 
given in the last column of the lower section of Table [3] 

From the calculations presented above, we conclude that 
only after the granulation effects have been properly taken into 
account, which implies reliably detecting the plateau in the line 
shift vs. equivalent width relation and accurately measuring the 
line shifts of the spectral lines that define the plateau, the uncer- 
tainty in the absolute radial velocity of a star is dominated by the 
err or in the gravitational redshift. 

iNidever et alj (|2002|) have determined the absolute radial ve- 
locities of a large sample of late-type stars, including HIP 96100, 
HIP 2331 1, and HIP 26779. Their approach combines the classi- 
cal method of cross-correlation with standard star templates with 
the very accurate relative radial velocity determination tech- 
niques used in spectroscopic exoplanet detection. Excluding sys- 
tematic errors, they claim to achieve a precision of about 20- 
30 m s _1 . For HIP 96100, HIP 23311, and HIP 26779, they 
derived 26.69, 21.55, and 1.21 km s _1 , respectively. The differ- 
ences with our values correspond to 0.44, 0.35 , and 0.67 km s " 1 , 
with our velocities being smaller in all cases. INidever et al.l es- 
timated their systematic error due to be about +0.30 km s _1 in 
K-dwarfs, which, if considered, would improve the agreement 
between our results and theirs, with the possible exception of 
HIP 26779, for which the large difference that remains suggests 
long-period variability. 

The high stability of the radial velocities of HIP 96100 and 
HIP 2331 1 make them good targets for radial velocity standard- 
ization. According to our calculations, the absolute radial ve- 
locities of these K-dwarfs are known wi th an accuracy of about 
60 m s -1 . Furthermore, for HIP 96100, lGrav et al.l d 19921) have 
shown that the line bisectors are time independent, which sug- 
gests that the line shifts (and therefore the granulation correc- 
tion to the absolute radial velocity) could also be constant, given 
that they are both observable manifestations of the same phe- 
nomenon. 

6. Summary and conclusions 

In this work, the signatures of granulation on the spectra of a 
small sample of K-dwarfs were detected using very high resolu- 
tion (R a 160, 000-210, 000), high signal-to-noise (S /N > 300) 
observations with a spectral coverage from 5580 to 7800 A. Not 
only were the data of very high quality, but also the reduction and 
post-reduction processing was carefully designed to preserve the 
relatively small effects of surface inhomogeneities. We were able 
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to determine relative radial velocity shifts with a mean accuracy 
of 12 m s _1 for pairs of spectra of the same object, without re- 
quiring the spectrograph to be placed in a vacuum chamber or su- 
perimposing iodine features on the stellar spectra. The accurate 
determination of relative shifts was crucial for coadding multi- 
ple exposures of the same object to increase the signal-to-noise 
ratio of the observations while keeping the data useful for the 
detection and measurement of the relatively weak granulation 
signatures, which require very high S IN. However, we cannot 
guarantee that these values correspond to the real relative veloc- 
ity variations of the stars because of shifts that could have been 
introduced by instrumental imperfections and/or environmental 
effects (e.g., variations of the air properties in the spectrograph 
room). In addition, we used a method of coaddition such that 
the interpolation distances in the resampling of frames before 
coadding was reduced by a factor of 2 compared to the standard 
method, thus reducing the impact of interpolation errors. 

Using our carefully processed spectra, we measured the bi- 
sectors of a large number of Fe i lines and averaged them out in 
groups of similar line-strength to minimize the impact of blends 
and noise. The observed mean bisectors show a characteristic 
C-shape, although for the strongest lines they have a shape that 
resembles only the lower part of the letter C (sometimes also 
referred to as a backslash shape: \), which spans up to about 
100 m s _1 , even though the exact value depends on the particu- 
lar K-dwarf star under consideration. We find that the stellar pro- 
jected rotational velocity and activity influence the shapes of the 
line bisectors, reducing their magnitudes compared to inactive 
stars with low V sin ; values. Although the effect of the projected 
rotational velocity is well understood, the details of the depen- 
dence on stellar activity, which we conclude it exists based on a 
correlation of level of chromospheric emission with line bisector 
span, remain open to investigation. 

Core wavelength shifts were also determined for all our K- 
dwarf sample stars and they were found to show the typical 
properties of stellar granulation; for example, the weakest lines, 
which are formed in deep regions of the photosphere, where the 
correlation between temperature and velocity fields is stronger, 
show the largest blueshifts (around -150 ms _1 ). Furthermore, 
the blueshifts decrease (i.e., the core wavelength shift, by defi- 
nition positive for redshift, increases) for stronger lines and, for 
some of our sample stars, they appear to reach a constant value 
for Fe i lines stronger than EW 100 mA. We interpret this 
plateau as line formation taking place in high photospheric lay- 
ers, where the correlation between intensity and velocity fields 
is very weak or null and therefore the convective line shifts are 
expected to be nearly zero. This behavior had been identified be- 
fore only in the solar spectrum. The detection of a plateau allows 
us to determine the zero point of the convective blueshifts, which 
is necessary to remove the large uncertainties due to granulation 
in the determination of absolute stellar radial velocities. In some 
of our K-dwarfs, the plateau could not be reliably detected, most 
likely due to extrinsic effects such as observational noise, scatter 
due to large rotational velocities, or stellar activity. 

Thanks to our quantitative understanding of the effects of 
granulation on the core wavelengths of Fe i lines in the spectra of 
K-dwarfs, we have been able to determine the absolute radial ve- 
locities of two stars that have non-variable radial velocities and 
a very well defined plateau in the convective blueshift vs. line 
strength relation (HIP 96100 and HIP 23311). These extremely 
accurate measurements make these stars, together with the Sun 
and asteroids with known orbits, the best radial velocity stan- 
dards available. In addition, we have measured radial velocities 
with similar precision for two stars that show long-term vari- 



ability: HIP 88601, which is the primary star of a visual binary 
system, and HIP 26779, a suspected radial velocity variable. 

The apparent non-detection of the plateau in other sample 
stars does not imply that granulation on their surfaces behaves 
differently, given that they suffer from other phenomena that pre- 
clude detection. Thus, in principle, it is possible to remove the 
errors due to granulation on measurements of absolute radial ve- 
locities of many K-type dwarf stars, which could be extremely 
useful for future massive spectroscopic surveys, where millions 
of K-dwarfs can be used to establish an accurate zero point for 
the stellar radial velocities. 
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Appendix A: Projected rotational velocities 

Most of our sample stars have one or more entries 
in The Catalog of Stel l ar Pr ojected Rotational Velocities 
dGlebocki & St awikowskil l2000h . which is a compilation of 
Vsin; values available in the literature and, therefore, deter- 
mined using a number of different techniques such as Fourier 
transform of line-profile (FTPL), cross-correlation (C-C), cali- 
brated line-width at half maximum (FWHM), and convolution 
with calculated rotational broadening (Conv). Table lA.fl lists the 
values given in this catalog for our sample of K-dwarfs (the orig- 
inal references are also provided). Entries giving only upper or 
lower limits as well as very old an d uncertain values have been 
excluded. Note that iFekell (11997) calibrated his method using 
measurements by D. F. Gray and his results ar e therefore n ot in- 
dependent from, for example, those given in iGravl (Il984al) . In 
addition, iFekelf s assumption of a constant macroturbulent ve- 
locity for K-dwarfs makes those results l ess re liable. The V sin/ 
values determined bv lValenti & Fischer! d2005l) are also given in 
Table \KJ\ 

By convolving line profiles predicted by a 3D K-dwarf 
model atmosphere with projected rotational velocity profiles, 
as it will be described in detail in Paper II, we determined a 
more accurate V sin i value for HIP 86400. The result is given in 
Table IA.1I Proper disk integrations were compared against the 
convolution method and the line profiles comput ed with both ap- 
proaches were found to agree remarkably well (Ramirez 2008, 
pp. 88-90). We did not perform similar calculations for the other 
sample stars because the line profiles are very sensitive to the 
stellar parameters, in particular effective temperature, and only 
HIP 86400 has parameters identical (within observational errors) 
to those of the 3D model. 



Appendix B: Chromospheric activity 

Many K-type dwarf stars reveal high levels of chromospheric 
activi ty in the form of e mission in the cores of very strong lines 
(e.g. jNoves eTaill 19841) . This may have an important impact on 
the shapes of spectral lines in general. In Fig. IB. II we show 
the spectral region around the Ca n H and K lines for all our 
sample stars except Pro cyon. The data were acqui red from the 
spectroscopic survey of Alle nde Prieto et ail (120041) but we have 
reduced the resolution to about R = 8, 000 for clarity in the il- 
lustration. The chromospheric activity index R'^ for each star is 
also given in Fig. IB. II This information is relevant to the discus- 
sion presented in Sect. |4] 
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Table A.l. Vsinz values reported in 


the literature for our K- 







dwarf sample stars. The value determined for HIP 86400 in 
Paper II of this work (TW) is also given. 



Appendix C: Orbital solution for the single-lined 
spectroscopic binary HIP 86400 

Comb ining our radial velocity data with those of iTokovininl 
d 1 99 lb . we were able to fit a spectroscopic orbit for the pri- 
mary star of the HIP 86400 system using the GaussFit package 
dJeffervs et al.l ll988). The resulting orbital parameters are given 
in Table IC. H and the model fit to the data is shown in Fig. IC.ll 

The mass function of the system is 0.00153 ± 0.00008M o , 
which implies a mass of about O.42M for the secondary, as- 
suming a mass of O.85M for the primary (determined from 
isochrones) and si n i = 0.32 (as measured astrometrically by 
iJancart et alj2005l) . This corresponds t o an M-dwarf star with an 
effective tempera ture around 3500 K (Baraffe & Chabrier 1996; 
Basri et aTll2000l) whose contribution to the flux is relatively 
small in the visible and near infrared (a secondary spectrum has 
not been detected yet). Using the spectral energy distributions 
calculated from Kurucz model atmospheres, we find that the sur- 
face flux from the M-dwarf in the HIP 86400 system is about 8% 
of the flux emitted by the K-dwarf in the spectral region of inter- 
est for this work. Considering that the flux received on Earth is 
proportional to the angular diameter of the star squared an d the 
radii of 0.4M o dwarfs are about 0.4# o dBeattv et al.ll2007l) . the 
contribution to the flux from the M-dwarf to our observed spec- 
trum of HIP 86400 is about (0.4/0.8) 2 X 8%= 2 %, given that the 
radius of the primary, the K-dwarf, is about O.87? . 




3920 3940 3960 3980 
Wavelength (A) 



3920 3940 3960 3980 
Wavelength (A) 



Fig. B.l. Medium resolution spectra of our K-dwarfs and the Sun 
in the region around the CanH and K line s. The activity index 
R' HY is shown in each panel, as given by ISaar & Brandenburg! 



; 1999 ), except for the cases marked with t he following sym bols: 
r *l =lRocha-Pinto & Maci 
= ISaar& Donahue! (1 19971) 



r*l = Rocha-Pinto & Maciel ( 19981) . [#] = lKing et al.l d2003l) . [$] 
9971). 



Parameter 



Value 



Units 



Orbital period (P) 


83.713 i 


.- 0.005 


days 


Time of periastron passage (T) 


2447723.54 ± 0.47 


JD 


Eccentricity (e) 


0.288 ± 


0.012 




Longitude of periastron (a>) 


135.99 


± 1.97 


degrees 


Semi-amplitude of RV curve (A*[) 


5.86 ± 


0.10 


km s" 1 


Velocity of the center of mass (y a ) 


22.095 i 


= 0.077 


km s -1 


Velocity of the center of mass (jt) 


21.194 d 


0.070 


km s" 1 


Mass function (f m ) 


0.00153 + 


0.00008 





Table C.l. Orbital parameters obtained by fitting the radial ve- 
locity data of HIP 86400. Two values for y are given; the first 
one, y a , corresponds to our observations, while yb was obtained 
for Tokovinin (1991) data. A systematic difference of about 
0.9 km s , with our velocities being larger, is present between 
the two data sets. The rest of the parameters were derived using 
the combined data set. The radial velocity curve of HIP 86400 is 
shown in Fig. IC.ll 
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Fig. C.l. Upper panel: radial velocities of HIP 86400 obtained 
from our observations. The dotted line is a model fit to the data. 
Lower panel: fit to our radial velocity data (black diamonds) and 
those of Tokovinin (1991; gray asterisks with error bars). The 
orbital parameters are given in Table IC.ll 



